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Abstract 


This  report  contains  a  description  and  performance  evaluation  of  the  matrix  calibration 
method  developed  by  the  author  in  1995  for  the  receive  antenna  array  in  the  high  frequency 
surface  wave  radar  (HFSWR)  at  Cape  Race,  Newfoundland.  The  evaluation  compares  this 
method  with  two  other  methods:  (a)  the  phase-only  calibration  method  being  used  and  (b)  an 
assumed  gain-and-phase  method  that  equalizes  the  complex  element  output  amplitudes  of 
signals  in  the  boresight  direction. 

The  results  showed  that  the  matrix  calibration  method  had  the  most  accurate  direction 
estimates,  the  narrowest  spectral  peaks,  the  best  resolution  of  closely  space  targets,  and  the 
lowest  signal-to-interference-plus-noise  ratio  (SINR)  threshold  for  target  detection.  The  phase- 
only  and  gain-and-phase  methods  had  false  target  tracks  and  direction-dependent  biases  in 
direction  estimates.  The  phase-only  method  also  had  a  line  splitting  problem. 

Based  on  the  difference  in  SINR  thresholds  for  target  detection,  the  matrix  calibration 
method  could  increase  the  present  maximum  nighttime  target  detection  range  of  the  HFSWR 
by  more  than  50%  in  situations  where  the  interference  is  highly  directional  and  extremely 
strong. 


Resume 


Ce  rapport  contient  une  description  et  une  evaluation  de  performance  de  la  methode  de 
calibration  de  matrice  developp6e  par  Tauteur  en  1995  pour  le  reseau  d’antenne  de  reception 
du  radar  decametriques  a  ondes  de  surface  (RDOS)  installee  a  Cap  Race,  Terre  Neuve.  La 
methode  de  calibration  de  matrice  est  comparee  avec  deux  autres  methodes  :  (a)  la  methode 
de  calibration  par  phase  seule  presentement  utilisee  dans  le  RDOS  et  (b)  une  methode  de 
phase  et  gain  qui  egalise  l’element  complexe  des  sorties  d’amplitude  des  signaux  dans  la 
direction  de  l’axe  d’antenne  de  pointage. 

Les  resultats  demontrent  que  la  methode  de  calibration  de  matrice  produit  les  mesures 
de  direction  les  plus  pr&ises,  les  cretes  spectrales  les  plus  etroites,  la  meilleure  resolution  de 
cibles  rapprochees,  et  le  seuil  de  detection  le  plus  bas  de  pour  le  rapport  signal  au  brouillage 
(SINR).  La  methode  de  phase  et  celle  de  phase  et  gain  produisent  des  fausses  traces  de 
cibles  et  des  erreurs  systematiques  dans  les  mesures  de  direction  dependant  de  la  direction  de 
la  cible.  La  methode  de  phase  a  aussi  des  problemes  de  division  des  lignes  spectrales. 

D’apres  la  difference  dans  le  seuil  de  detection  pour  le  SINR,  la  methode  de 
calibration  de  matrice  peut  augmenter  la  portee  maximum  de  detection  de  nuit  du  RDOS  par 
plus  de  50%  pour  des  situations  ou  le  brouillage  est  tres  directionnel  et  extremement  fort. 
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Executive  summary 


The  Defence  Research  Establishment  Ottawa  (DREO)  and  Raytheon  Canada  Limited 
have  jointly  built  a  high  frequency  surface  wave  radar  (HFSWR)  at  Cape  Race, 
Newfoundland,  for  coastal  surveillance  on  a  24-hours/7-days  basis.  This  radar  is  located 
behind  a  cliff  facing  the  Atlantic  Ocean  and  is  designed  to  detect  and  track  targets,  e.g., 
surface  ships  and  low-flying  aircraft,  below  the  microwave  radar  horizon.  It  can  work  in 
isolation  or  as  part  of  a  larger  surveillance  system  that  incorporates  other  sensors  such  as 
satellites,  line-of-sight  radars,  and  other  HFSWRs. 

The  HFSWR  has  been  operational  since  summer  1998.  A  problem  that  has  not  yet 
been  solved  is  the  detection  of  long-range  targets  at  night,  when  extremely  strong  skywave 
interference  signals  are  present.  These  signals  originate  from  distance  sources.  They  can  reach 
the  HFSWR,  because  the  signal-absorbing  D  layer  of  the  ionosphere  is  absent  at  night.  During 
daylight  hours,  this  radar  has  no  difficulties  detecting  large  ships  at  ranges  up  to  450  km  from 
the  radar.  At  night,  the  maximum  range  drops  dramatically  to  less  than  200  km.  Many  attempts 
to  increase  this  range  had  been  made.  None  could  extend  the  maximum  range  beyond  200  km. 

The  skywave  interference  was  one  of  several  problems  identified  in  1994  when  a 
prototype  HFSWR  was  operating  at  Cape  Race.  A  matrix  calibration  method  for  the  receive 
antenna  array  of  the  radar  was  developed  in  1995  by  this  author  to  mitigate  their  detrimental 
effects.  This  method  effectively  uses  sixteen  sets  of  complex  calibration  coefficients  at  once  to 
compensate  for  the  unequal  gain  and  phase  responses  at  the  array  elements,  mismatches  in 
cables  and  electronics  between  the  elements  and  the  receivers,  small  errors  in  element  locations, 
mutual  coupling  among  the  elements,  imperfect  ground  screens,  irregular  shoreline,  signal 
diffraction  at  the  cliff,  and  interaction  between  the  array  and  other  local  objects.  The  reasoning 
was  that  detection  of  targets  outside  the  interference  direction  was  still  possible,  if  the  array  was 
accurately  calibrated  to  reduce  the  leakage  of  interference  power  outside  the  interference 
direction. 

This  report  presents  a  description  and  a  performance  evaluation  of  the  matrix 
calibration  method,  using  a  set  of  calibration  data  collected  in  April  2000.  The  evaluation 
compares  the  matrix  method  with  two  other  methods:  (a)  the  single-set  phase-only  calibration 
method  being  used  and  (b)  an  assumed  single-set  gain-and-phase  method  that  equalizes  the 
complex  element  output  amplitudes  of  signals  in  the  boresight  direction.  Some  notable 
observations  on  targets  within  55°  of  the  boresight  direction  are  listed  below. 

(a)  The  leakage  of  signal  power  outside  the  source  direction  were  0.04%  and  1.29%  in  data 
sections  with  high  and  low  signal-to-interference-plus-noise  ratios  (SINRs),  respectively,  in 
calculations  with  the  matrix  calibration  method.  They  were  between  5.20  and  11.7%  in  the 
single-set  phase-only  method,  and  between  0.21  and  5.73%  in  the  single-set  gain-and- 
phase  method. 

(b)  The  errors  in  source  direction  estimates  were  usually  the  smallest  in  the  matrix  method  and 
the  largest  in  the  phase-only  method. 
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(c)  The  errors  in  source  direction  estimates  had  a  component  linearly  dependent  on  the  source 
azimuth  in  the  phase-only  and  gain-and-phase  methods. 

(d)  The  signal  peaks  in  the  directional  spectra  of  the  array  snapshots  were  the  narrowest  in 
calculations  with  the  matrix  method.  They  were  significantly  broader  in  the  phase-only 
method.  Those  in  the  gain-and-phase  method  were  just  as  narrow  as  the  peaks  in  the  matrix 
method  for  signals  in  the  boresight  direction.  They  increased  to  about  half  the  widths  of  the 
phase-only  method  as  the  source  moved  towards  the  ±50°  directions. 

(e)  The  power  levels  of  false  peaks  were  less  than  -40  dB  of  the  true  signal  in  the  matrix 
method.  They  were  above  -30  dB  and  frequently  over  -20  dB  in  the  phase-only  method. 
Those  in  the  gain-and-phase  method  are  usually  between  -30  and  -20  dB. 

(f)  The  false  peaks  did  not  form  false  target  tracks  in  the  matrix  method.  They  frequently 
formed  false  tracks  in  the  phase-only  and  gain-and  phase  methods. 

(h)  Line  splitting  (two  direction  estimates  separated  by  less  than  one  beamwidth  in  high- 
resolution  direction-finding,  one  signal  source  actually  present)  was  very  common  in 
calculations  with  the  phase-only  method. 

(i)  The  SINR  threshold  in  detecting  targets  outside  the  region  of  extremely  strong  interference 
was  the  lowest  in  the  matrix  method.  It  was  approximately  17  dB  higher  in  the  phase-only 
method  and  5  to  8  dB  higher  in  the  gain-and-phase  method. 

Based  on  the  approximate  1 7  dB  difference  in  SINR  threshold  for  target  detection,  the  matrix 
calibration  method  could  increase  the  present  maximum  nighttime  target  detection  range  of 
the  HFSWR  by  more  than  50%  in  situations  where  the  interference  is  highly  directional  and 
extremely  strong. 


Hung,  EKL.  2001.  "A  calibration  procedure  for  an  HFSWR  receive  array".  DREO  TR 
2001-103.  Defence  Research  Establishment  Ottawa. 
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Sommaire 


Le  Centre  de  Recherche  de  la  Defense  d’Ottawa  et  Raytheon  Canada  Limited  ont 
construit  conjointement  un  radar  decimetrique  a  onde  de  surface  (RDOS)  a  Cap  Race,  Terre 
Neuve.  Ce  radar  a  ete  construit  pour  la  surveillance  cotiere  24  heures  par  jour  et  7  jours  par 
semaine.  Le  RDOS  est  situee  sur  une  falaise  en  face  de  1’  ocean  Atlantique  et  est  con9u  pour 
la  detection  et  la  poursuite  des  cibles  se  deplaqant  en  dessous  de  1‘horizon  radar  microonde  tel 
que  des  navires  et  des  avions  volant  a  basse  altitude.  Le  RDOS  peut  operer  seule  ou  etre 
integrer  dans  un  systeme  de  surveillance  plus  grand  qui  incorpore  d’autres  recepteurs  tel  que 
les  satellites,  les  radars  de  lignes  visuelles,  et  d’autres  RDOS. 

Le  RDOS  a  ete  operationnel  depuis  l’ete  1998.  La  detection  des  cibles  k  longue 
portee  durant  la  nuit,  quand  le  brouillage  ionospherique  est  intense,  est  un  probleme  qui  n’a 
pas  ete  resolu  juste  qu’a  maintenant.  Ce  brouillage  ionospherique  provient  de  sources 
distantes.  Ces  signaux  lointains  de  brouillages  peuvent  atteindre  le  RDOS  durant  la  nuit  car  la 
couche  ionospherique  D  qui  absorbe  les  signaux  d6cametrique  est  absente  durant  la  nuit. 

Durant  le  jour,  le  RDOS  n’a  pas  de  difficulte  a  detecter  de  grands  navires  jusqu’a  450km  de 
portee.  Durant  la  nuit,  la  portee  maximum  diminue  dramatiquement  a  seulement  200km. 
Plusieurs  essais  ont  ete  faits  pour  augmenter  cette  portee,  mais  aucun  n’a  ete  capable  de 
l’augmenter  au-dela  de  200km. 

Le  brouillage  ionospherique  est  un  des  problemes  identifies  en  1994  quand  un 
prototype  de  RDOS  operait  a  Cap  Race.  Une  methode  de  calibration  de  matrice  pour 
l’antenne  receptrice  du  radar  a  ete  developpee  en  1995  par  l’auteur  pour  reduire  leurs  effets 
negatifs.  Cette  methode  utilise  16  ensemble  de  coefficients  de  calibration  complexe  pour 
compenser  le  gain  et  les  reponses  inegal  des  elements  du  reseau  d’antenne,  le  facteur 
d’inegalite  des  cables  et  de  l’electronique  entre  les  elements,  les  petites  erreurs  dans  la 
position  des  elements,  le  couplage  mutuelle  entre  les  elements,  les  imperfections  de  l’ecran 
de  terre,  l’irregularite  de  la  ligne  de  rivage,  la  diffraction  du  signal  sur  la  falaise,  et 
1’ interaction  entre  le  reseau  et  les  objets  locales.  Le  raisonnement  est  que  la  detection  des 
cibles  a  l’exterieur  du  brouillage  est  encore  possible,  si  le  reseau  est  precisement  calibrer  pour 
reduire  la  puissance  de  brouillage  a  l’exterieur  de  la  direction  de  brouillage. 

Ce  rapport  contient  une  description  et  une  evaluation  de  performance  de  la  methode 
de  calibration  de  matrice  utilisant  un  ensemble  de  calibration  collecter  en  avril  2000.  La 
methode  de  calibration  de  matrice  est  compare  avec  deux  autres  methodes  :  (a)  la  methode 
de  calibration  par  phase  seule  presentement  utilisee  dans  le  RDOS  et  (b)  une  methode  de 
phase  et  gain  qui  egalise  1’ element  complexe  des  sorties  d’ amplitude  des  signaux  dans  la 
direction  de  l’axe  d’antenne  de  pointage.  Quelques  observations  notables  sur  les  cibles 
situees  a  moins  de  55  degres  de  1’axe  de  pointage  de  l’antenne  sont  listees  ci  dessous. 

(a)  La  puissance  du  signal  a  l’exterieur  de  la  direction  de  la  source  etait  de  0.04%  et 
1.29%  dans  la  section  des  donnees  ayant  des  rapports  de  signal  k  brouillage  (SINR) 
faible  et  fort,  respectivement,  dans  les  calculs  utilisant  la  methode  de  calibration  de 
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matrice.  Elies  etaient  entre  5.2  et  1 1.7%  dans  la  methode  de  phase  seul,  et  entre 
0.21  et  5.73  dans  la  methode  de  calibration  de  gain  et  phase. 

(b)  Les  erreurs  dans  les  mesures  de  direction  etaient  habituellement  les  plus  petites 
avec  la  methode  de  calibration  de  matrice  et  etaient  les  plus  grandes  dans  la 
methode  de  phase  seule. 

(c)  Les  erreurs  dans  les  mesures  de  la  direction  de  la  source  ont  une  composante 
lineaire  qui  depend  de  l’azimut  pour  les  methodes  de  calibration  de  phase  seule  et 
aussi  de  phase  et  gain. 

(d)  Les  cretes  des  signaux  dans  le  spectre  directionnelles  des  donnees  instantanees  sont 
les  plus  etroites  quand  on  utilise  la  methode  de  calibration  de  matrice.  Les  cretes 
sont  beaucoup  plus  larges  quand  on  utilise  la  methode  de  calibration  de  phase 
seule.  Les  cretes  obtenues  avec  la  methode  de  phase  et  gain  sont  aussi  etroites  que 
celles  obtenues  avec  la  methode  de  calibration  de  matrice.  Cependant  les  cretes 
elargissent  jusqu’a  la  moitie  de  la  largeur  de  celles  obtenues  avec  la  methode  de 
calibration  de  phase  seule  lorsque  la  source  se  deplace  loin  de  1’axe  de  ‘antenne 
vers  55  degres. 

(e)  Les  niveaux  de  puissance  sont  moins  que  -40dB  comparee  au  signal  reel  avec  la 
methode  de  calibration  de  matrice.  Ils  sont  au-dessus  de  -30dB  et  souvent  au- 
dessus  de  -20db  dans  la  methode  de  calibration  de  phase  seule.  Ceux  de  la 
m&hode  de  calibration  de  gain  et  phase  sont  habituellement  entre  -30  et  -20dB. 

(f)  Les  fausses  cretes  n’ont  pas  produit  de  fausses  traces  de  poursuite  avec  la  methode 
de  matrice.  Cependant,  ils  ont  souvent  formee  des  fausses  traces  de  poursuite  avec 
la  methode  de  phase  seul  et  aussi  la  methode  de  calibration  de  phase  et  gain. 

(g)  La  division  des  lignes  spectrales  (deux  mesures  de  directions  separees  par  moins 
que  une  largeur  de  faisceau  d’antenne  dans  la  radiogoniometrie  haute  resolution,  un 
signal  source  actuellement  present)  sont  tres  commun  avec  la  methode  de 
calibration  de  phase  seule. 

(h)  Le  seuil  de  SINR  pour  la  detection  de  cibles  a  l’exterieur  de  la  region  de  brouillage 
ionospherique  intense  est  le  plus  bas  avec  la  methode  de  calibration  de  matrice.  II 
est  approximativement  1 7dB  plus  grand  avec  la  methode  de  phase  seule  et  de  5dB 
a  8dB  avec  la  methode  de  phase  et  gain. 

D’apres  la  difference  de  17dB  dans  le  seuil  de  detection  pour  le  SINR,  la  methode  de 
calibration  de  matrice  peut  augmenter  la  portee  maximum  de  detection  de  nuit  du  RDOS  par 
plus  de  50%  pour  des  situations  oil  le  brouillage  est  tres  directionnel  et  extremement  fort. 


Hung,  EKL.  2001.  "A  calibration  procedure  for  an  HFSWR  receive  array".  DREO  TR 
2001-103.  Defence  Research  Establishment  Ottawa. 
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1 .  Introduction 


The  Defence  Research  Establishment  Ottawa  (DREO)  and  Raytheon  Canada  Limited 
have  jointly  developed  a  high  frequency  surface  wave  radar  (HFSWR)  at  Cape  Race, 
Newfoundland,  for  coastal  surveillance  on  a  24-hours/7-days  basis.  This  radar  is  located 
behind  a  cliff  feeing  the  Atlantic  Ocean.  It  has  a  uniform,  linear  receive  array  with  sixteen 
monopole  doublets  spaced  at  a  distance  of  33.33  meters.  The  radar  is  designed  to  detect  and 
track  targets,  e.g.,  surface  ships  and  low-flying  aircraft,  below  the  microwave  radar  horizon.  It 
can  work  in  isolation  or  as  part  of  a  larger  surveillance  system  that  incorporates  other  sensors, 
including  satellites,  line-of-sight  radars,  and  other  HFSWRs. 

Detection  of  long-range  targets  at  night  requires  detection  of  weak  target  reflectors  in 
the  presence  of  strong  skywave  interference.  The  signals  in  this  interference  originate  from 
distant  sources.  They  can  travel  thousands  of  kilometres  to  reach  the  HFSWR  at  night, 
because  the  signal-absorbing  D  layer  of  the  ionosphere  is  only  present  during  the  day.  To 
separate  the  target  from  interference,  it  is  necessary  to  use  high-resolution  direction-finding 
methods  such  as  the  Capon  minimum  variance  method  [1-3]  and  the  least  squares 
autoregressive  method  (LSAR)  [4,5],  or  interference  suppression  methods  such  as  the  sample 
matrix  inverse  method  [6,7]  and  the  adaptive  sidelobe  cancellation  method  [8],  or  both.  These 
methods  require  accurately  calibrated  antenna  arrays  [9-14]. 

The  skywave  interference,  together  with  a  few  others,  was  identified  by  Dr.  Tim 
Coyne,  in  1994.  He  concluded  that  a  proper  calibration  of  the  receive  antenna  array  was 
essential  to  mitigate  their  detrimental  effects.  The  reasoning  was  that  detection  of  targets 
outside  the  interference  direction  is  still  possible,  if  the  interference  is  highly  directional  and  if 
the  array  is  accurately  calibrated  to  reduce  the  leakage  of  interference  power  outside  the 
interference  direction.  A  matrix  calibration  method  for  the  HFSWR  was  developed  by  this 
author  in  1995.  It  effectively  used  sixteen  sets  of  complex  coefficients  at  once  to  compensate 
for  the  unequal  gain  and  phase  responses  at  the  array  elements,  mismatches  in  cables  and 
electronics  between  the  elements  and  the  receivers,  small  errors  in  element  locations,  mutual 
coupling  among  the  elements,  imperfect  ground  screens,  irregular  shoreline,  signal  diffraction  at 
the  cliff,  and  interaction  between  the  array  and  other  local  objects.  Because  the  HFSWR  was  not 
operational  at  the  time,  the  method  was  evaluated  with  data  from  four  other  sources.  All  results 
indicated  significant  improvements  in  target  detection  after  calibration  [15-18]. 


The  HFSWR  at  Cape  Race  has  been  fully  operational  since  summer  1998.  However,  it 
is  still  being  calibrated  with  a  single-set  phase-only  method  that  only  addresses  the  unequal 
phase  responses  at  the  array  elements.  During  daylight  hours,  this  radar  has  no  difficulties 
detecting  large  ships  up  to  450  km  from  the  radar.  At  night,  the  maximum  range  drops 
dramatically  to  less  than  200  km.  Many  attempts  to  increase  this  range  had  been  made.  All  used 
the  phase-only  calibration  method.  In  agreement  with  what  Dr.  Coyne  had  concluded,  none 
could  extend  the  maximum  range  beyond  200  km. 
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This  report  presents  a  description  and  a  performance  evaluation  of  the  matrix 
calibration  method,  using  a  set  of  calibration  data  collected  in  April  2000.  The  evaluation 
compares  the  matrix  method  with  two  other  methods:  (a)  the  single-set  phase-only  calibration 
method  being  used  and  (b)  an  assumed  single-set  gain-and-phase  method  that  equalizes  the 
complex  element  output  amplitudes  of  signals  in  the  boresight  direction.  In  the  results,  the 
matrix  calibration  method  has  the  most  accurate  direction  estimates,  the  narrowest  spectral 
peaks,  the  best  resolution  of  closely  space  targets,  and  the  lowest  signal-to-interference-plus- 
noise  ratio  (SINR)  threshold  for  target  detection;  the  phase-only  and  gain-and-phase  methods 
show  false  target  tracks  and  direction-dependent  biases  in  direction  estimates;  and  the  phase- 
only  method  also  has  a  line  splitting  problem. 

Based  on  the  difference  in  SINR  thresholds  for  target  detection,  the  matrix  calibration 
method  could  increase  the  present  maximum  nighttime  target  detection  range  of  the  HFSWR 
by  more  than  50%  in  situations  where  the  interference  is  highly  directional  and  extremely 
strong. 


In  the  following,  Section  2  describes  the  matrix  calibration  method  developed  in 
1995.  Section  3  discusses  the  calibration  data  and  the  calibration  matrix.  Section  4  presents 
the  results  of  performance  evaluations  and  Section  5  contains  the  conclusions.  There  are  two 
appendices.  Appendix  A  describes  the  computation  of  signal  direction  estimates  with  the 
LSAR  method  and  Appendix  B  describes  the  computation  of  directional  spectra  with  the 
Capon  method. 
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2.  Calibration  method 


A  brief  description  of  the  calibration  method  developed  in  1995  is  presented  here.  It 
includes  the  signal  model  used,  the  conditions  that  must  be  satisfied  and  the  computation  of 
the  calibration  matrix. 


2.1  Signal  model 

The  receive  antenna  array  has  N  elements.  A  pilot  source  is  present  and  the  source 
azimuth  direction  is  0m.  An  ideal  snapshot  of  the  outputs  of  the  array  elements  is  given  by 


Xm  OCmJlm , 

where  a™  is  the  complex  signal  amplitude  and  am  is  the  array  steering  vector  for  the  source 
direction.  This  expression  assumes  a  perfectly  calibrated  array,  a  pure  tone  signal,  no 
interaction  among  the  array  elements,  no  noise,  and  no  interference.  In  the  convention  used, 
angle  0m  is  measured  clockwise  from  the  array  boresight  direction. 

The  corresponding  measured  array  snapshot  is  given  by  ym,  with 

ym  0Xm  +  e,n  . 

Matrix  B  is  NxN.  It  depends  on  the  gain  and  phase  responses  at  the  array  elements,  mismatches 
in  cables  and  electronics  between  the  elements  and  the  receivers,  small  errors  in  element 
locations,  mutual  coupling  among  the  elements,  imperfect  ground  screens,  irregular  shoreline, 
signal  diffraction  at  the  cliff,  and  interaction  between  the  array  and  other  local  objects.  Vector 
e,,,  denotes  the  sum  of  interference  and  noise. 

Let  {y,,y2,...,yM}  be  a  set  of  M  array  snapshots  measured  with  source  directions 
{0i, 02, —,0m},  respectively.  The  calibration  matrix  calculated  with  this  set  is  identified  here  as 
the  NxN  matrix  C  that,  together  with  the  estimates  of  {x!,x2,...,Xm},  minimises  the  objective 
function 


S  = 


M 

Z  |Cym“Xmf 

m—  l _ 

M 

m— i 


subject  to  the  uniqueness  condition 


C„  =  1.0 . 


(3) 


(4) 
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Vector  Xm  is  an  estimate  of  xm  and  will  be  calculated  as  the  projection  of  Cym  on  am.  The 
denotes  the  estimate  and  Cpq  denotes  the  pq-th  coefficient  of  C. 

In  later  discussions,  the  superscripts  T  and  H  denote  the  transpose  and  conjugate- 
transpose,  respectively;  the  set  {yim,y2m,-,  yNm},  for  example,  denotes  the  N  components  of 
ym;  and  y™  defined  as 

ycm  =  Cym ,  (5) 

denotes  the  calibrated  snapshot  calculated  with  ym. 

It  may  be  noted  that  C  has  the  singular  value  decomposition 


C  =  UXVH,  (6) 

where  U=(Ui,u2,...,Un)  and  V=(v,,v2,...,vN)  are  matrices  with  orthonormal  column  vectors,  and 
Z  is  a  diagonal  matrix  with  singular  values  {c2,,o22,...,Cn}.  Substitution  into  (5)  gives, 

N 

ycm=  X1o2„(vHnym)un .  (7) 

n— 1 

This  expression  shows  that  C  effectively  calculates  the  projection  of  ym  on  each  vn,  i.e.,  the 
projection  (v”„ym)vn,  and  then  replaces  this  projection  by  (v*iym)Un-  Therefore,  using  Eq.  (5)  to 
calculate  yin  is  equivalent  to  using  N  sets  of  calibration  coefficients  at  once  to  calibrate  the 
array  snapshot  ym.  Each  set  has  N  complex  amplitudes. 

The  computation  of  C  is  a  pure  mathematical  procedure.  The  function  S  depends  on 
the  signal  amplitude  estimates  cXi  to  aM  through  x,  to  Xm,  respectively.  In  the  special  case  of 
zero  noise  and  interference,  the  tm  in  each  ym  is  a  null  vector,  CB=(3In,  Cym=Xm,  xm=(3xm,  and 
S=0.  The  P  is  a  complex  scalar  and  IN  is  an  NxN  identify  matrix. 


2.2  Conditions 

The  following  conditions  are  assumed  to  be  satisfied. 

1 .  The  array  elements  are  narrowband; 

2.  The  process  that  converts  an  ideal  snapshot  to  a  measured  snapshot  can  be  approximated 
by  a  matrix; 

3.  The  calibration  data  are  generated  by  moving  a  pilot  source  over  the  range  of  azimuth 
directions  covered  by  the  HFSWR; 

4.  The  number  of  array  snapshots  in  the  calibration  data  is  sufficiently  large;  and 

5.  A  preliminary  calibration  has  been  carried  out,  so  that  one  can  obtain  rough  estimates  of 
the  complex  signal  amplitudes  from  the  array  snapshots. 

Conditions  3  and  4  are  needed  for  tolerance  to  unknown  interference,  including 

diffuse  and  specula  multipath.  Together,  they  attempt  to  destroy  the  correlation  between  the 
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pilot  signal  and  any  interference  present  in  the  environment  by  effectively  randomising  their 
phase  difference  in  the  data.  They  also  cause  each  fixed-direction  interference  source  to 
behave  like  spatially  distributed  interference  by  forcing  the  angular  separation  between  it  and 
the  pilot  signal  to  spread  over  a  wide  range  of  values.  The  number  of  snapshots  required  in 
Condition  4  is  the  number  required  for  C  to  reach  the  asymptotic  limit  as  this  number 
increases. 


2.3  Calculation  of  calibration  matrix 

The  calculation  of  C  is  based  on  the  following  observations. 
A.  If  Cm  is  an  estimate  of  C  and 


Zm  Ci-lYm » 

a  least  square  estimate  of  the  signal  amplitude  is 


*  _  amZn, 
CMn  * 


famf  ’ 

and  an  estimate  of  xm  is 


Xm  OlntSlm  . 


(8) 

(9) 


(10) 


B.  Given  the  NxM  matrices  Xand  Y,  where  X=(xi,x2,...,Xm)  and  Y-(yi,y2,...,yM),  one  can 
get  a  new  estimate  of  C  by  solving  the  equation 


CjY  =  X 


(11) 


as 

Ci  =  (XXh)(YXh)-'  ,  (12) 

and  imposing  the  uniqueness  condition  (4)  on  Ci.  The  inverse  (YXH)''  exists,  because 
of  Conditions  3  and  4. 


C.  Given  snapshot  ym  and  the  array  steering  vector  an,  for  source  direction  0™,  the 
initialisation 


Co=diag{ami/ymi,am2/ym2,...,amN/ymN} ,  (13) 

where  diag{...}  denotes  a  diagonal  matrix,  is  equivalent  to  a  crude  calibration  that 

equalises  the  element  output  gains  and  phases  for  signals  with  source  direction  0m. 
The  set  {ami,aII12,...,ainN}  denotes  the  components  of  am. 
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D. 


Let 


Dj  =  Cj  -  Cm 


(14) 


and 


A;  =  max{Ppq|:  p,q=l,2,...,N} .  (1:>) 

One  can  terminate  the  iteration  if  (a)  Aj<£,  where  £  is  a  small  value,  or  (b)  i-'max- 
where  imax  is  the  maximum  number  of  iterations  specified. 

There  are  six  steps  in  the  procedure.  The  Sm  in  Step  3  is  the  value  of  S  in  the  (i-l)th 
iteration. 


Step  1 . 

(a)  Specify  e. 

(b)  Specify  imax* 

(c)  Let  i=l. 

(d)  Calculate  C0  with  Eq.  (13)  and  the  snapshot  ym  with  0m  closest  to  0°. 


Step  2. 

Calculate  {zuz2,  ...,Zm}  with  Eq.  (8). 

Calculate  {ii,x2,  ...,Xm}  with  Eqs.  (9)  and  (10). 


Step  3. 

Calculate  Sm  with  Eq.  (3),  {ii,x2,  —  ,Xm}>  and  {Z|,Z2,...,Zm}- 
Calculate  Cj  with  Eq.  (12). 


Step  4. 


Impose  the  uniqueness  condition  given  by  Eq.  (4)  on  Cj,  i.e.,  the  condition  (Cj)it  1  • 


Step  5. 

Calculate  Ai  with  Eqs.  (14)  and  (15). 


Step  6. 

If  Aj<e  or  i>imax. 

Terminate  the  iterations  and  identify  Cj  as  C. 

Else, 

Replace  i  by  i+1  and  return  to  Step  2. 

In  Step  3,  vector  z™  is  the  product  C,.iym  needed  to  calculate  SM.  The  calculation  of 
Sj.,  can  be  omitted  if  the  sole  objective  is  to  check  if  S,  converges.  From  Eq.  (3),  S,  converges 
if  Q  converges. 
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3.  Experiment  and  calibration  matrix 


This  section  discusses  the  calibration  experiment,  data  preprocessing,  the  calibration 
matrix,  and  effective  calibration  coefficients  for  given  source  directions. 

The  calibration  matrix  calculated  is  only  useful  for  calibrating  array  snapshots  in  which 
all  signal  and  interference  directions  are  more  than  one  beamwidth  from  the  ends  of  the  azimuth 
range  -80°  to  58.7°,  i.e.,  for  signals  and  interference  direction  inside  the  range  -54.0°  to  42.7°. 


3.1  Experiment 

The  receive  antenna  array  in  the  HFSWR  at  Cape  Race  is  located  behind  a  23  m  high 
cliff  facing  the  Atlantic  Ocean.  It  has  a  uniform  linear  receive  antenna  array  with  sixteen 
monopole  doublet  elements,  i.e.,  N=16,  spaced  at  a  distance  of  33.33  m.  Figure  1  is  a  map  of  the 
radar  site.  It  shows  the  location  of  the  antenna  array,  a  hut  that  houses  the  radar  receivers  and 
other  hardware,  a  lighthouse  northeast  of  the  array,  and  the  shoreline  near  the  lighthouse.  The 
coastline  is  irregular.  The  land  behind  the  cliff  was  boggy  and  had  many  streams  and  ponds. 

A  pilot  source  on  board  a  44  foot  fishing  boat  was  used  to  generated  the  calibration 
signals  at  the  receive  array.  Originally,  the  boat  was  to  sail  from  the  0=90°  endfire  direction  of 
the  antenna  array  to  the  0=-9O°  endfire  direction  at  a  range  of  28  km  from  the  array.  Due  to 
poor  weather  conditions,  a  section  of  the  track  southwest  of  Cape  Race  was  omitted  and  the 
boat  actually  sailed  from  0=58.7°  to  0=-9O°. 

The  pilot  source  consisted  of  an  ICOM  IC-735  HF  transceiver  free  running  at  3.225 
MHz  and  transmitting  at  a  nominal  power  level  of  50  watts.  This  transceiver  was  powered  by 
an  ICOM  power  supply,  which  in  turn  was  powered  by  the  boat's  110-volt  generator.  It  fed  a 
BIRD  power  coupler  and  meter,  to  provide  a  measure  of  the  forward/reflected  power,  and 
then  an  antenna  matching  unit,  model  MFJ  901B,  from  MFJ  Inc.  The  unit  was  connected  to  a 
22  foot  HF  whip  antenna  mounted  on  the  boat's  rigging.  The  lower  end  of  the  whip  was 
approximately  25  feet  above  the  waterline  and  the  grounds  of  all  components  were  attached  to 
the  boat's  internal  ground  system. 

At  the  antenna  site,  the  output  of  each  array  element  was  connected  by  cable  to  a 
receiver  and  heterodyned  to  25  kHz.  Quadrature  sampled  snapshots  of  the  array  element  outputs 
were  taken  in  sets,  identified  as  pulses,  of  64  snapshots  per  pulse  and  at  a  rate  of  62.5  pulses  per 
second.  Groups  of  eight  pulses  were  then  summed  to  produce  7.8125  preaveraged  pulses  per 
second.  The  64  snapshots  in  each  pulse  were  sampled  at  10  microsecond  intervals,  or  100  kHz. 

Data  collection  started  at  approximately  10:00  Newfoundland  standard  time  and 
ended  at  16:00.  There  were  instances  when  the  boat  was  asked  to  "heave  to"  and  standby  in 
position,  so  that  downloading  of  collected  data  could  be  carried  out  at  the  hut.  A  DGPS 
(differential  global  positioning  system)  provided  the  time  and  boat  position.  These  values 
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were  logged  by  training  a  video  camera  on  the  DGPS  display.  As  a  precaution,  the  values 
were  also  hand  logged  every  5  minutes. 


3.2  Data  preprocessing 

The  calibration  data  collected  had  186077  preaveraged  pulses  with  64  array  snapshots 
per  pulse,  16  complex  amplitudes  per  snapshot,  and  8  bytes  per  complex  amplitude. 
Preprocessing  was  necessary,  because 

(a)  The  pulses  were  unevenly  spaced  in  source  azimuth; 

(b)  Each  receiver  had  a  zero  frequency  bias  amplitude  that  depended  on  the  receiver; 

(c)  There  were  many  strong  interfering  spikes  in  the  data;  and 

(d)  The  file  had  over  1 .5  gigabytes  of  data  and  had  to  be  reduced  to  a  manageable  size. 

The  error  amplitudes  in  (b)  originated  from  a  flaw  in  receiver  design  and  the  spikes  in  (c)  were 
generated  by  an  ionosonde  that  was  programmed  to  turn  on  for  several  minutes  just  before  every 
GMT  hour.  The  ionosonde  was  turned  off  at  the  controls  in  the  hut  during  the  calibration 
experiment.  Unknown  to  the  operator,  it  was  still  sounding  the  ionosphere  at  hourly  intervals. 

The  data  preprocessing  procedure  was  as  follows. 

Step  1.  Bandpass  filtering 

The  outputs  of  each  element  in  each  preaveraged  pulse  were  transformed  with  a  64- 
point  fast  Fourier  transform  (FFT)  to  produce  a  range  FFT  spectrum  with  64  complex 
amplitudes.  The  amplitude  in  the  25  kHz  signal  bin  was  extracted  from  each  element  to 
produce  a  bandpass  filtered  snapshot. 

Step  2.  Doppler  filtering 

The  intermediate  file  constructed  with  the  bandpass  filtered  snapshots  was  divided  into 
overlapping  groups  with  a  window  256  snapshots  wide.  The  window  moved  down  the 
file  in  steps  of  32  snapshots.  Each  group  of  256  snapshots  was  then  converted  to  a 
Doppler  filtered  snapshot  with  a  256-point  FFT. 

Step  3.  Removal  of  snapshots  with  large  errors  in  source  direction  estimates 

The  Doppler  filtered  snapshots  were  calibrated  with  the  phase-only  calibration 
coefficients  and  the  source  direction  estimate  0  in  each  calibrated  snapshot  was  then 
calculated  with  the  LSAR  method  [4,5].  All  snapshots  with  |0m-8n,t>6o  were  deleted. 

Step  4.  Thinning 

The  snapshots  at  the  output  of  Step  3  were  selectively  deleted  to  make  the  remaining 
snapshots  approximately  equally  spaced  in  azimuth. 

The  preprocessed  calibration  data  file  was  denoted  by  {yi,y2,—,yM}-  It  had  3470 
snapshots,  i.e.,  M=3470,  spaced  approximately  0.032°  apart  in  the  range  -8O.O°<0m<58.7Oo.  For 
reasons  to  be  given  later,  all  snapshots  with  6m<-80.0°  were  deleted.  The  computation  of  source 
direction  estimate  is  given  in  Appendix  A. 
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Step  1  was  designed  to  remove  the  zero  frequency  bias  at  the  receiver  outputs  and  to 
increase  the  SINR  at  the  array  elements.  The  25  kHz  signal  bin  was  bin  17  in  the  64-point  FFT, 
because  the  array  snapshots  in  each  pulse  were  sampled  at  100  kHz.  In  Step  2,  the  256  complex 
bandpass  filtered  amplitudes  of  each  element  in  the  group  were  converted  to  a  Doppler  spectrum 
with  a  256-point  FFT.  Next,  the  Doppler  bin  with  the  largest  power  in  the  spectrum  of  the  first 
element  was  identified  as  the  signal  bin  in  the  Doppler  spectra  of  all  elements.  Then,  the 
complex  amplitude  in  this  bin  was  extracted  from  the  spectrum  of  each  element  and  used  to 
construct  a  Doppler  filtered  snapshot.  Step  3  was  designed  to  remove  Doppler  filtered  snapshots 
with  unusually  low  SINRs.  Before  this  step,  there  were  265  snapshots  with  6m<-80.0°.  After  this 
step,  only  61  remained.  Because  over  75%  were  removed,  the  remaining  snapshots  had 
questionable  reliabilities  and  were  removed  as  well.  No  snapshots  with  -75.0°<em<58.70o  were 
removed  and  very  few  with  -8O.O°<0m<-75.O0  had  |Mml>60.  Step  4  was  necessary,  because  the 
array  snapshots  were  very  unevenly  spaced  in  azimuth. 

The  reason  for  the  large  number  of  snapshots  with  0m-6ml>6o  for  0m<-8O.O°  could  be  the 
large  land  mass  northeast  of  the  receive  antenna  array.  It  could  also  be  strong  mutual  coupling 
among  the  array  elements  for  these  source  directions.  Indications  of  mutual  coupling  can  be 
found  in  Figures  7b  and  7c  in  the  performance  evaluation  to  be  presented  later. 

Figure  2  shows  a  typical  set  of  N,  where  N=16,  64-point  FFT  power  spectra  calculated 
with  the  outputs  of  the  array  elements  in  one  preaveraged  pulse.  To  avoid  unnecessary 
complications  from  the  ionosonde  spikes,  the  nominal  source  azimuth  direction  was  chosen  to 
be  near  0.0°,  so  that  it  was  more  than  6.0°  from  the  nearest  spike.  There  are  two  peaks  in  the 
spectrum  of  each  element.  The  peak  at  25  kHz  is  the  signal  peak  and  the  one  at  0  kHz  peak  is 
the  bias  at  the  receiver.  Table  1  shows  the  bias  and  signal  powers  at  the  elements.  These  powers 
are  between  1245  and  4102  units  at  the  25  kHz  signal  bins  and  between  0.17  and  19.76  umts  at 
the  0  kHz  bin.  The  conversion  factor  from  power  units  to  signal  field  strengths  at  the  array 
elements  was  unavailable. 

Figure  3  shows  a  typical  256-point  Doppler  power  spectrum  calculated  in  Step  2  with 
the  complex  signal  amplitudes  in  the  25  kHz  signal  bin  of  the  64-point  FFT  spectra  in  Step  1 .  In 
this  example,  the  signal  amplitudes  were  calculated  with  the  outputs  of  Element  8  and  the 
nominal  source  direction  was  also  near  0.0°.  The  spectrum  has  a  band  with  many  peaks  that  are 
15  dB  or  more  above  the  noise  background.  The  bandwidth  is  larger  than  1.0  Hz.  A  study  of  the 
data  showed  that  the  shape  of  the  band  changed  with  time  and  that  the  band  position  drifted 
slowly  in  frequency.  Figure  4  shows  the  256-point  Doppler  power  spectrum  obtained  when  the 
complex  bias  amplitudes  in  the  0  kHz  bin  were  used  instead.  It  has  a  strong  peak  at  0.0  Hz. 

The  bandwidth  of  more  than  1.0  Hz  in  Figure  3  was  unexpectedly  large,  because  a 
monotone  source  was  specified  for  the  experiment.  The  drift  in  band  position  was  expected, 
because  the  pilot  source  on  board  the  fishing  boat  was  free  running  and,  therefore,  not 
synchronized  with  the  oscillator  in  the  receivers. 

Figure  5a  shows  the  average  power  at  the  array  elements  in  the  256  preaveraged  pulses 
used  in  Steps  1  and  2  to  generate  one  preprocessed  snapshot.  There  are  strong  spikes  at 
0  =48.2°,  27.9°,  9.4°,  -18.9°  and  -48.8°.  From  the  records  on  time  and  boat  position,  these  are 
interference  spikes  generated  by  the  ionosonde,  because  they  are  one  hour  apart.  The  time  of  the 
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rightmost  spike  is  GMT  13:59:46  (NST  10:29:46)  and  the  time  of  the  leftmost  spike  is  GMT 
1 7:59:46  (NST  14:29:46).  The  peaks  are  not  equally  space  in  direction,  because  the  boat  did  not 
move  at  a  constant  angular  speed  along  the  track.  Figure  5b  shows  the  average  power  at  the 
elements  in  the  preprocessed  snapshots.  The  preprocessing  in  Steps  1  and  2  had  removed  the 
spikes.  It  had  also  reduced  the  average  power  outside  the  spike  regions  by  approximately  10 
units. 


3.3  Calibration  matrix 

The  array  snapshots  in  the  preprocessed  calibration  data  file  were  used  to  calculate  C, 
using  the  procedure  at  the  end  of  Section  2.3.  In  Step  1,  the  assignments  were  £=0.00001  and 
imax=50.  In  Step2,  the  array  steering  vector  in  Eq.  (9)  was  calculated  as 

am  =  (l,expU2jtdsin0m],...,  exp[j2jid(N-l)sin0m])T ,  (16) 

where  d=0.35552=(33.33/93.75)  was  the  element  spacing  in  wavelengths.  In  this  expression, 
the  first  array  element  was  the  element  farthest  to  the  north.  In  Step  3,  the  matrix  product  YXH 
in  Eq.  (12)  was  highly  illconditioned,  because  there  were  no  snapshots  with  9m  outside  the 
range  -80.00°  to  58.70°.  Consequently,  the  inverse  (YX”)'1  was  replaced  by  the  Moore- 
Penrose  pseudo-inverse  calculated  with  the  singular  value  decomposition  method.  Only  the 
thirteen  strongest  eigenvectors  were  retained,  because  this  number  produced  the  best  results  in 
the  evaluation  of  the  calibration  method.  In  the  intermediate  results,  the  Aj  in  Eq.  (15)  dropped 
quickly  from  1.20016  at  i=l  to  0.02197  at  i=2.  It  then  decreased  gradually  to  0.00056  at  i=50. 
The  corresponding  value  of  Sj.i  dropped  from  0.95132  at  i=l  to  0.00712  at  i=2  and  to  0.00677 
at  i=7.  It  then  increased  very  slowly  to  0.00689  at  i=29  before  decreasing  again  to  0.00682  at 
i=50.  Table  2  shows  the  magnitudes  of  the  coefficients  {CM}  in  the  upper  left  8x7  comer  of 
C.  Evidently,  C  was  not  symmetric  and  the  magnitudes  in  the  main  diagonal  were  different.  A 
check  on  the  validity  of  Condition  4  in  Section  2.2  could  not  be  made,  because  the 
rank(C)=13<N. 

Matrix  C  was  effectively  the  asymptotic  limit  of  Q  for  extremely  large  in**.  This 
property  was  verified  by  calculating  a  second  C,  denoted  by  C,  using  £=0.00000  and 
W=200.  Matrix  C  looked  identical  to  C  when  the  first  three  decimal  places  of  the  matrix 
coefficients  were  compared.  It  also  produced  near  identical  results  in  calculations.  Therefore, 
for  all  practical  purposes,  C  and  C  were  the  same  and  C  had  effectively  converged  to  the 
asymptotic  limit.  A  rigorous  proof  of  convergence  is  beyond  the  scope  of  this  report. 

Matrix  C  was  not  unique.  The  C's  calculated  with  different  C0's  were  not  expected  to 
produce  the  same  results  from  snapshots  with  source  directions  less  than  one  beamwidth  from 
the  ends  of  the  data  azimuth  range  -80.00°  to  58.70°.  However,  all  were  expected  to  produce 
essentially  the  same  results  if  the  directions  were  more  than  one  beamwidth  from  the  ends. 
These  properties  were  verified  by  calculating  a  C'  with  a  Co  derived  from  the  phase-only 
calibration  coefficients  used  to  calibrate  the  HFSWR.  In  agreement  with  the  above  remarks, 
matrices  C  and  C  looked  different  but  produced  near  identical  directional  spectra  from 
snapshots  with  0m  more  than  one  beamwidth  from  the  ends.  The  C0  in  the  calculation  of  C 
was  the  phase-only  calibration  matrix  CP  in  the  performance  evaluation  in  Section  4. 
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Matrix  C  was  not  unique,  because  YXH  was  highly  illconditioned  and  the  weak 
eigenvectors  of  YXH  depended  on  C0  as  well  as  the  ym’s  with  0m's  less  than  one  beamwidth 
from  the  ends.  These  weak  eigenvectors  determined  the  performance  of  C  for  snapshots  with 
source  directions  up  to  one  beamwidth  inside  the  data  azimuth  range  or  totally  outside  it.  All 
C's  were  expected  to  produce  essentially  the  same  results  for  source  directions  more  than  one 
beamwidth  from  the  ends,  because  the  calibration  data  determined  the  strong  eigenvectors  of 
YXh.  These  eigenvectors  in  turn  determined  the  performance  of  C  for  snapshots  with  these 
source  directions. 

Matrix  C  was  only  useful  for  calibrating  array  snapshots  in  which  all  signal  and 
interference  directions  were  more  than  one  beamwidth  from  the  ends  of  the  azimuth  range  -80° 
to  58.7°.  Because  azimuth  direction  measured  in  beam  widths  from  boresight  was  given  by 

0b  =  Ndsin0 ,  (1?) 

C  was  only  useful  for  calibration  array  snapshots  with  -54.O°<0m<42.70. 


3.4  Effective  calibration  coefficients 

The  calibrated  snapshot  ycm  and  the  uncalibrated  snapshot  ym  are  related  by  Eq.  (5), 
i  e.,  ycm=Cym.  It  follows  from  this  relation  that  the  effective  calibration  coefficients  for  the 
components  of  ym  are  given  by 


Y«>  = 


n=l,2,...,N , 


(18) 


where  {yWc2m,...,yCNm}  and  {y^y^.-.jNm}  are  the  components  of  ycm  and  y™  respectively. 
These  coefficients  are  used  to  multiply  the  outputs  of  the  array  elements.  They  depend  on  0m 
through  ym.  Table  3  contains  the  magnitudes  and  phases  of  effective  calibration  coefficients 
(Yi,Y2,...,Yn}  for  0m=O°  and  0m=-6O°.  In  the  same  table  are  the  phase-only  coefficients  used  to 
calibrate  the  HFSWR.  For  the  sake  of  convenience  in  comparison,  all  three  sets  of  coefficients 
have  been  normalised  to  give  71=1.0,  so  that  abs(Yi)=l  and  angle(yi)=0.  There  are  several 
notable  observations.  Firstly,  the  coefficients  calculated  with  Eq.  (18)  are  highly  dependent  on 
0m.  In  the  special  case  of  the  eighth  array  element,  i.e.  n=8,  for  example,  the  magnitude  is  0.74 
for  0m=O°  and  0.57  for  0m=-6O°.  Secondly,  for  a  given  0,,,,  the  magnitudes  and  phases  of  the 
coefficients  depend  on  the  array  elements.  Thirdly,  the  magnitudes  and  phases  of  the  phase- 
only  coefficients  are  different  from  those  calculated  with  Eq.  (18). 
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Element 

Bias  Power 

Signal  Power 

(0  kHz) 

{25  kHz) 

1 

2.17 

1245 

2 

1.98 

1257 

3 

0.98 

2668 

4 

1.24 

2668 

5 

0.17 

2785 

6 

2.78 

1942 

7 

5.55 

3011 

8 

19.76 

2469 

9 

0.63 

4102 

10 

3.35 

2503 

11 

1.04 

3456 

12 

0.47 

2985 

13 

3.80 

2349 

14 

2.74 

3010 

15 

4.48 

1991 

16 

1.53 

1487 

Table  1:  Bias  power  and  signal  power  in  the  range  FFT  spectra  of  array  elements  in  a  preaveraged 
pulse  with  source  azimuth  dose  to  the  0.0°  boresight  direction.  The  source  frequency  was  3.225  MHz. 
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p 

q=i 

m 

3 

4 

5 

6 

7 

1 

1.000 

0.057 

0.099 

0.068 

0.084 

0.009 

0.028 

2 

0.073 

0.833 

0.133 

0.140 

0.181 

0.052 

0.046 

3 

0.074 

0.104 

0.601 

0.154 

0.090 

0.127 

0.032 

4 

0.017 

0.102 

0.120 

0.500 

0.179 

0.109 

0.039 

5 

0.055 

0.073 

0.136 

0.150 

0.504 

0.180 

0.114 

6 

0.033 

0.071 

0.092 

0.111 

0.158 

0.513 

0.138 

7 

0.021 

0.038 

0.022 

0.050 

0.074 

0.091 

0.474 

8 

0.055 

0.068 

0.072 

0.012 

0.063 

0.167 

0.181 

Table  2:  Magnitudes  of  some  coefficients  {Cpq}  at  the  top  left  comer  of  calibration  matrix  C. 
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Magnitude 

Phase 

n 

P 

11 

0 

0 

O 

O 

1 

II 

CD 

P 

0m=o° 

em=-60° 

P 

1 

1.00 

1.00 

1.00 

0 

0 

0 

2 

1.11 

0.89 

1.00 

-161 

-162 

-146 

3 

0.86 

0.49 

1.00 

44 

41 

52 

4 

0.71 

0.51 

1.00 

-92 

-114 

-100 

5 

0.71 

0.57 

1.00 

79 

66 

99 

6 

0.85 

0.60 

1.00 

-63 

-64 

-48 

7 

0.66 

0.42 

1.00 

134 

131 

151 

8 

0.74 

0.57 

1.00 

-15 

-30 

-12 

9 

0.58 

0.45 

1.00 

170 

159 

-162 

10 

0.72 

0.51 

1.00 

-29 

-35 

-5 

11 

0.64 

0.44 

1.00 

122 

119 

134 

12 

0.68 

0.43 

1.00 

-99 

-101 

-91 

13 

0.71 

0.50 

1.00 

100 

90 

107 

14 

0.68 

0.46 

1.00 

-108 

-114 

-104 

15 

0.90 

0.56 

1.00 

29 

23 

27 

16 

0.90 

0.52 

1.00 

151 

145 

151 

Table  3:  Magnitudes  and  phases  of  some  calibration  coefficients  {Y1.Y2 . Yn}-  Columns  identified  with 

6m=0°  and  -60°  were  calculated  with  Eq.  (17).  Those  identified  with  P  are  used  in  phase-only  calibration. 
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Figure  1:  Site  map  of  the  high  frequency  surface  wave  radar  at  Cape  Race. 
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Figure  2:  A  set  of  sixteen  64-point  FFT  power  spectra  calculated  with  the  outputs  of  the  array 
elements  in  one  preaveraged  pulse  measured  with  source  azimuth  direction  dose  to  0.0°. 
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Figures  5a  and  5b:  Average  power  at  array  elements  in  the  256  preaveraged  pulses  used  to  produce 
one  preprocessed  snapshot  (Figure  5a,  top)  and  average  power  at  array  elements  in  the 
preprocessed  snapshots  (Figure  5b,  bottom). 


4.  Performance  evaluation 


Matrix  C  and  two  other  matrices,  denoted  by  Cp  and  Cop,  were  used  in  Eq.  (5)  to 
calibrate  the  preprocessed  snapshots  {yi,y2,-»yM}-  The  calibrated  snapshots  were  then  used  to 
study  the  leakage  of  signal  power  outside  the  signal  direction,  errors  in  source  direction 
estimates,  directional  power  spectra,  signal  direction  estimates,  target  resolution,  and  target 
detection  in  the  presence  of  strong  interference.  Cp  and  Cgp  were  diagonal  matrices  calculated 
with  the  coefficients  identified  with  P  and  9m=0°,  respectively,  in  Table  3,  e.g.. 


CP  =  diag{Yi,Y2,— ,Yn}  ■  (19) 

Using  CP  is  equivalent  to  using  the  single-set  phase-only  calibration  method  being  used  to 
calibrate  the  HFSWR  and  using  Cgp  is  equivalent  to  using  the  assumed  single-set  gain-and- 
phase  method  that  equalizes  the  complex  element  output  amplitudes  of  signals  in  the  boresight 
direction. 

Due  to  limitations  in  the  calibration  data,  special  emphases  were  given  to  results  in 
the  azimuth  range  -75°  to  55°.  Whenever  there  was  a  choice  of  signal  source  directions,  these 
directions  were  restricted  to  the  range  -54.0°  to  42.7°,  because  C  was  most  effective  for 
signals  in  this  range. 


4.1  Leakage  power  outside  signal  direction 

From  Eqs.  (1),  (2)  and  (5),  one  gets 
y  m  =  Cym 
=  C(amBam  +  em) 

=  Pn>(am+bm)  +  e’m ,  (20) 

where  pm(am+bm)  and  e’m  represent  the  products  c^CBa,*  and  Ce™,  respectively.  In  this 
expression,  pm  is  a  complex  scalar  representing  the  calibrated  signal  amplitude  at  the  array 
elements,  and  bm  is  a  vector  generated  by  calibration  error.  By  definition,  bm  is  orthogonal  to 
aIT„  i.e. ,  ambm  bmarTi  0. 

In  Eq.  (20),  the  component  of  e’m  parallel  to  the  signal  direction  is  given  by 

the  component  of  e'm  orthogonal  to  the  signal  direction  is 
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Cm_|_  ©m  •  Cm//  ? 


(22) 


the  component  of  ycm  orthogonal  to  the  signal  direction  is 


yCmX  =  Pmbm  +  e'mi.  ,  (23) 

and  the  fraction  of  total  power  in  ycm  outside  the  signal  direction  and  averaged  over  all  array 
elements  is 


L(9)  =  •  <24> 

The  signal  term  in  the  calibrated  snapshot  ycm  is  Pm(am+bm)  and  the  term  generated  by 
calibration  error  is  pjbm.  When  the  SINR  and  the  error  (5mbm  are  both  so  large  that  pm(am+bm)[ 
»|e'in|  and  iMU^jeUt  one  gets  yCm=Pm(am+bm),  yCmi=M>m  and 


LIM. 

SINR,  pb|— >oo 


L(6)  = 


|bmf 

K+brnf 


(25) 


This  limit  is  smaller  if  the  antenna  array  is  more  accurately  calibrated,  because  a  more 
accurate  calibration  method  produces  a  smaller  Ibmt2.  It  usually  has  a  weak  dependence  on  the 
source  direction,  because  bm  depends  on  the  source  direction  through  the  <XmCBam  term  in  Eq. 
(20). 


The  curve  C  in  Figure  6  shows  the  values  of  L(0)  calculated  with  C  and  every  tenth 
snapshot  in  the  preprocessed  calibration  data  file.  It  has  spikes  close  to  0m=48.2o,  27.9°,  9.4°, 
and  -18.9°.  These  spike  directions  are  also  the  directions  of  the  first  four  ionosonde  spikes  in 
Figure  5a.  Outside  the  spike  regions,  the  values  are  between  -34.0  and  -18.9  dB, 
corresponding  to  0.04  and  1.29%,  respectively,  of  total  power.  There  are  no  signs  of 
ptnbmf»|e'mj.l)  because  the  largest  and  smallest  values  of  L(0)  differ  by  more  than  15  dB. 
Besides,  from  Figure  5b,  the  regions  of  lower  and  larger  values  of  L(0)  also  correspond  to  the 
regions  of  larger  and  lower  average  powers,  respectively,  at  the  array  elements.  Curve  Cp  was 
calculated  with  CP.  Here,  the  values  of  L(0)  vary  between  -12.9  dB  (5.20%)  and  -9.3  dB 
(1 1 .7%),  thus  giving  a  difference  of  only  3.6  dB  between  the  two  limits.  This  small  difference 
of  3.6  dB  is  a  clear  indication  of  |(5mbm|»|e'mil  Curve  Cgp  was  calculated  with  CGP.  The 
values  of  L(0)  are  between  -26.8  dB  (0.21%)  and  -12.4  dB  (5.73%).  With  the  exception  of  a 
small  region  close  to  0m=O°,  they  are  larger  than  those  calculated  with  C  by  more  than  5  dB. 

The  leakage  azimuth  spectrum  of  a  pilot  source  is  also  the  leakage  spectrum  of  the 
signals  in  nighttime  interference.  These  interfering  signals  are  usually  over  30  dB  stronger 
than  the  long-range  target  signals  the  HFSWR  is  designed  to  detect.  Figure  6  indicates  serious 
target  detection  problems  when  CP  is  used  to  calibrate  the  array  snapshots.  Suppose  a  target 
signal  is  -30  dB  relative  to  the  interference;  the  leakage  from  the  interference  is  equal  to  the 
lower  value  of  5.2%;  and  this  leakage  is  uniformly  distributed  over  all  directions.  Here,  the 
average  leakage  power  is  52  times  the  target  power  at  the  array  elements.  Coherent  integration 
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can  increase  the  SINR  by  a  factor  16,  which  is  the  number  of  elements  in  the  receive  array. 
Even  with  this  integration  gain,  the  average  leakage  power  in  the  target  direction  is  3.25 
(=52/16)  times  the  target  power.  Therefore,  the  probability  of  detecting  the  target  is  very  low. 
Usually,  the  leakage  produces  one  or  more  false  target  peaks  at  unpredictable  azimuth 
directions.  The  probability  of  detecting  the  true  target  in  the  presence  of  these  false  targets  is 
even  lower. 


4.2  Errors  in  direction  estimates 

The  LSAR  method  was  used  to  calculate  a  set  of  five  candidate  signal  direction 
estimates  from  every  tenth  array  snapshot  with  -75°<0m<550.  The  estimate  with  the  largest 
power  was  identified  as  the  source  direction  estimate  in  the  snapshot.  Given  an  estimate  0m,  the 
error  in  this  estimate  was  calculated  as 


Err(em)  =  0m-em.  (26) 

The  calculation  of  candidate  signal  direction  estimates  and  the  powers  associated  with  the 
estimates  is  given  in  Appendix  A. 

Figure  7a  shows  the  errors  obtained  when  C  was  used  to  calibrate  the  snapshots.  The 
errors  are  plotted  in  dot  mode.  A  horizontal  line  has  been  added  to  study  their  deviations  from 
zero.  In  the  range  -6O°<0m<55°,  most  errors  have  magnitudes  less  than  0.2°.  Figures  7b, 
calculated  with  Cp,  generally  has  errors  with  much  larger  magnitudes.  These  errors  have  a 
component  linearly  dependent  on  the  source  azimuth.  They  also  increase  rapidly  in  magnitude 
as  0m  falls  below  -70°.  Figure  7c,  calculated  with  Cqp,  has  errors  that  are  generally  between 
those  in  Figures  7a  and  7b  in  magnitudes.  Like  those  in  Figure  7b,  the  errors  have  a  linear 
component  and  the  error  magnitude  increase  rapidly  as  0m  falls  below  -70°. 

The  linear  components  of  errors  in  Figures  7b  and  7c  are  consistent  with  the  presence  of 
weak  mutual  coupling  among  the  array  elements  for  source  directions  between  -60°  and  55°. 
The  rapid  increase  for  0m<-7O°  could  be  caused  by  a  rapid  increase  in  coupling. 


4.3  Directional  power  spectra  of  calibration  data 

The  Capon  minimum  variance  method  [1-3]  was  used  to  calculate  the  directional 
power  spectra  of  every  tenth  snapshot  with  -75o<0m<55o.  These  spectra  were  defined  as 

p.Q.  _  Power  arriving  at  array  elements  from  spectral  direction  0  ^7) 

'  '  Total  power  at  array  elements 

A  brief  description  of  the  calculation  of  Capon  directional  spectra  is  given  in  Appendix  B. 

Figure  8a  shows  the  dependence  of  P(0)  on  source  azimuth  0m  when  C  was  used  to 
calibrate  the  preprocessed  snapshots.  It  has  a  strong  source  track  from  the  lower  left  comer, 
where  0=0m=-75°,  to  the  upper  right,  where  0=0m=55°.  For  each  6m,  there  is  at  least  one 
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spectral  azimuth  region  with  P(0)<-4O  dB.  The  figure  also  indicates  generally  larger  values  of 
P(9)  for  0<-5O°  and  0>5O°.  Figure  8b  is  calculated  with  CP.  The  source  track  is  significantly 
broader;  the  values  of  P(0)  outside  the  main  beam  are  significantly  higher;  and  there  are  large 
interference  patches  up  to  the  -17  dB  level  all  over  the  place.  Figure  8c  is  calculated  with  C<a>. 
The  source  track  is  narrowest  in  the  0m=O°  azimuth  region.  Outside  this  region,  the  width 
increases  as  the  deviation  of  0m  from  0°  increases.  At  0m=±55°,  the  widths  are  approximately 
half  of  those  calculated  with  CP.  The  figure  also  has  interference  patches  up  to  the  -17  dB 
level.  Compared  with  Figure  8b,  the  patches  are  fewer  in  number  and  are  farther  away  from 
the  source  azimuth. 

Figures  9a  to  9c  each  contain  the  spectra  of  six  snapshots  with  0m=4O°  in  Figures  8a  to 
8c,  respectively.  In  Figure  9a,  P(0)  falls  rapidly  to  approximately  -36  dB  as  the  deviation  of  0 
from  0m  increases  to  about  10°.  Left  of  the  signal  peak,  P(0)  continues  to  fall  and  is  below  -40 
dB  in  the  region  -5O°<0<25O.  For  0<-5O°,  P(0)  is  generally  larger  than  -40  dB  and  some  values 
are  close  to  -36  dB.  In  Figure  9b,  the  signal  peaks  are  significantly  broader  than  those  in  Figure 
9a  and  the  values  of  P(0)  outside  the  signal  region  are  up  to  20  dB  higher.  There  are  also  false 
peaks  in  the  region  0=-7°.  This  region  is  inside  an  interference  patch  in  Figure  8b.  In  Figure  9c, 
the  widths  of  the  signal  peaks  and  the  values  of  P(0)  outside  the  signal  region  are  between 
those  in  Figures  9a  and  9b.  There  are  false  peaks  in  the  region  0=2°. 

The  slightly  higher  background  for  0<-5O°  in  Figures  8a  and  9a  probably  originated 
from  diffuse  signal  scattering  at  the  cliff,  which  is  northeast  of  the  receive  array  in  Figure  1. 

The  directional  spectrum  of  a  pilot  source  at  azimuth  0m  is  also  the  directional 
spectrum  of  a  strong  sky  wave  interference  at  interference  azimuth  0i,  with  0i=0m-  Patches  of 
interference  up  to  the  -17  dB  level  in  Figure  8b  are  signs  of  total  disaster  in  target  detection  at 
night  when  CP  is  used.  These  patches  generate  false  targets  and  bury  the  long-range  targets  at 
once. 


4.4  Direction  estimates  calculated  with  calibration  data 

The  LSAR  method  was  used  to  calculate  five  candidate  signal  direction  estimates 
from  every  tenth  preprocessed  snapshot.  The  four  with  the  strongest  powers  were  then  plotted 
in  dot  mode  and  the  fifth  was  discarded  to  avoid  overcrowding.  Figures  10a,  10b  and  10c 
show  the  results  obtained  by  calibrating  the  snapshots  with  C,  CP  and  CgP,  respectively. 
Figure  10a  gives  the  correct  source  track,  which  is  a  straight-line  joining  the  lower  left  comer, 
where  0=0m=-75°,  with  the  upper  right,  where  0=0m=55°.  Figures  10b  and  10c  give  source 
tracks  that  are  slightly  tilted.  The  tracks  are  also  slightly  bent  at  the  lower  left  comer.  Both 
figures  have  false  tracks.  Compared  with  Figure  10b,  the  false  tracks  in  Figure  10c  are  fewer 
in  number  and  are  farther  away  from  the  source  track. 


The  shapes  of  source  tracks  in  Figures  10b  and  10c  agree  with  the  errors  in  source 
direction  estimates  in  Figure  7b  and  7c.  The  false  tracks  are  inside  the  interference  patches  in 
Figure  8b  and  8c.  The  figures  do  not  contain  the  same  number  of  dots,  because  some  direction 
estimates  are  non-physical  and  some  are  outside  the  range  in  the  figures. 
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4.5  Target  resolution 

Target  resolution  was  studied  with  a  synthetic  data  file  containing  a  moving  target 
and  a  fixed  target  with  the  same  strength.  The  snapshots  in  the  file  were  denoted  by  the  set 
{Zi,z2,...,Zm}  and  were  generated  as 

Zn,  =  exp[<|>m]ym  +  exp[<t>mF]yF ,  m=l,2,...,M  ,  (28) 


where 


ym  = 


ym 

tyJ 


(29) 


was  preprocessed  snapshot  ym  normalised  to  unit  magnitude,  yF  was  the  ym  with  source 
direction  0m=0F,  0f  was  the  direction  of  the  fixed  target,  and  <)>„,  and  <tw  were  random  phases. 


A  set  of  526  synthetic  snapshots  was  generated  with  0f=3O°  and  snapshots  with  0m 
between  -50°  and  -10°.  Figures  11a,  lib  and  11c  contain  the  Capon  directional  spectra 
obtained  by  using  C,  CP  and  Cgp,  respectively,  to  calibrate  the  synthetic  snapshots.  The  two 
target  tracks  are  the  narrowest  in  Figure  1  la  and  the  broadest  in  Figure  11c.  These  tracks  are 
resolved  in  Figure  1  la  if  pm-0p|  >5°.In  contrast,  they  are  not  resolved  in  Figures  1  lb  and  1  lc 
unless  0m-©F|>lOo.  There  is  a  false  target  track  between  0=-5°  and  0=15°  in  Figure  1  lc. 

Some  Capon  directional  spectra  with  0m-0F  are  extracted  from  Figures  1  la  to  1  lc  and 
plotted  in  Figures  12a  to  12c,  respectively.  Each  spectrum  has  a  strong  peak  close  to  0F.  The 
peaks  in  Figure  12a  are  the  narrowest  and  those  in  Figure  12b  are  the  broadest.  Figure  12c  has 
false  peaks  with  P(0)>-23  dB  in  the  region  0*10°.  These  false  peaks  are  inside  the  false  track 
in  Figure  1  lc. 

Figures  13a  to  13c  show  the  two  strongest  candidate  signal  direction  estimates 
calculated  with  the  calibrated  snapshots  and  the  LSAR  method.  In  agreement  with  the  results 
in  Figure  1  la  and  1  lc,  the  targets  are  resolved  in  Figure  13a  and  13c  if  Pm-fM  ^5°  and  pm-0F| 
>10°,  respectively.  However,  in  total  disagreement  with  the  results  in  Figure  1  lb,  the  targets 
in  Figure  13b  are  resolved  for  all  values  of  0m,  including  the  special  case  of  0m=0F.  There  is  a 
false  target  in  Figure  13c. 

The  contradicting  results  on  target  resolution  in  Figure  lib  and  13b  indicate  that 
something  catastrophic  had  happened.  The  LSAR  method  attempts  to  fit  an  array  snapshot 
with  sinusoids,  with  each  sinusoid  associated  with  a  different  signal  direction.  When  a 
snapshot  has  a  directional  spectrum  with  a  broad  and  strong  peak,  the  method  may  fit  the 
snapshot  with  more  than  one  sinusoids  close  in  directions.  Line  splitting  occurs  when  it 
generates  two  direction  estimates  less  than  one  one-way  array  beamwidth  apart  and  only  one 
signal  source  is  actually  present.  One  one-way  array  beamwidth  in  array  steering  direction  0, 
with  0=0f,  is  approximately  10.3°  (=  50.8°  /  Ndcos0)  for  0f=3O°. 

Broadening  of  a  strong  signal  peak  occurs  if  there  are  errors  in  antenna  calibration.  In 
the  case  of  a  poorly  calibrated  antenna  array,  the  peak  could  be  so  broad  that  the  LSAR 
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method  generates  more  than  one  direction  estimates  close  to  the  true  source  direction.  Line 
splitting  had  occurred  when  CP  was  used  to  calibrate  the  array  snapshots  and  the  LSAR 
method  then  used  to  calculate  signal  direction  estimates.  The  case  of  0m=0F  corresponded  to  a 
single  strong  signal  source  present.  One  should  not  have  obtained  two  estimates  close  to  0p. 

Whenever  there  is  a  possibility  of  line  splitting,  the  existence  of  two  estimates  close 
in  direction  does  not  automatically  imply  the  existence  of  two  targets  close  in  direction.  All 
claims  of  targets  close  in  directions  must  therefore  be  verified  with  information  from  other 
sources,  e.g.,  prior  knowledge  of  two  targets  present. 

Other  examples  of  line  splitting  can  be  found  in  Figure  10b  in  regions  with  source 
azimuth  0m--35°,  -5°  and  42°. 


4.6  Target  detection  in  presence  of  strong  interference 


Target  detection  in  the  presence  of  strong  interference  was  studied  with  a  synthetic 
file  containing  K=401  snapshots,  a  target  at  azimuth  direction  0r=-3O°,  and  J=6  interference 
signals  with  interference  directions  close  to  40°.  The  snapshots  were  generated  as 

zk=pkexp[<])k]yT 

+  exp[<>ki]yii  +  exp[<}>k2]yi2  +  ...  +  exp[<}>u]yu ,  k=l,2,...,K,  (30) 

Parameter  pk  was  a  real  scalar  related  to  the  target  SINR  in  by 


SINR  =  -  (31) 

The  set  {4b$u4ia-4u}  were  random  phases,  yT  was  the  normalised  array  snapshot  ym  with 
0m=8r,  and  the  set  {yn,yi2,...,yu}  was  a  set  of  normalised  snapshots  with  source  directions 
close  to  40°.  The  directional  spectra  of  these  J  interference  snapshots  are  given  in  Figures  9a 
to  9c. 


Two  sets  of  results  are  presented.  The  first  set  is  calculated  with  the  LSAR  method 
and  the  second  is  calculated  with  the  discrete  Fourier  transform  (DFT)  method.  The  results  of 
the  DFT  method  are  included,  because  this  method  is  being  used  in  target  detection  with  the 
HFSWR  at  night. 

Figure  14a  shows  the  Capon  directional  spectra  obtained  by  using  C  to  calibrate  the 
synthetic  snapshots.  The  tracks  at  0=-3O°  and  40°  are  generated  by  the  target  and  interference 
respectively  The  power,  width  and  azimuth  fluctuations  in  these  tracks  reflect  the  use  of 
random  phases  on  the  right  hand  side  of  Eq.  (30).  As  SINR  increases  from  -50  dB,  faint  target 
peaks  begin  to  appear  at  SINR=-42  dB.  The  peaks  become  clearly  visible  for  SINR>-35  dB. 
Figure  14b  shows  the  spectra  calculated  with  Cp.  The  target  and  interference  peaks  are  wider. 
With  very  few  exceptions,  the  target  peak  is  not  visible  unless  SINR>-23  dB.  There  are  false 
target  tracks  at  0=-55°  and  -5°.  Figure  14c  was  calculated  with  CqP.  Here,  the  peak  widths  are 
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between  those  calculated  with  C  and  Cp;  the  target  peaks  are  faint  at  S1NR  -36  dB  and  are 
clearly  visible  for  SINR>-30  dB;  and  there  are  signs  of  a  false  track  at  0=5  . 

Figures  15a,  15b  and  15c  each  contains  ten  Capon  directional  spectra  in  Figures  14a, 
14b  and  14c,  respectively.  The  SINR’s  were  =-40,  -23  and  -33  dB,  respectively,  and  were 
chosen  to  make  50%  of  the  spectra  in  each  figure  possess  a  peak  close  to  the  target  direction. 
Besides  the  five  target  peaks.  Figure  15a  has  four  broad  peaks  with  P(0)<-32  dB  at  0=-6O°.  It 
also  has  a  spectrum  with  a  significantly  stronger  signal  peak  and  a  very  distorted  interference 
peak.  Figure  15b  has  a  false  peak  with  P(0)=-19  dB  at  0=-8  in  nine  of  the  ten  spectra.  There 
are  signs  of  broad  peaks  at  0=-6O°.  All  values  of  P(0)  for  0<-5O°  exceed  -22  dB.  Figure  1 5c 
has  a  false  peak  with  P(0)=-24  dB  at  0=5°  in  every  spectrum.  The  spectra  are  essentially  flat 
for  0<-5O°. 

The  peaks  at  0=-6O°  in  Figure  15a  were  probably  generated  by  diffuse  scattering  of 
signals  at  the  cliff  northeast  of  the  receive  antenna  array.  These  scattered  signals  were  too 
weak  to  produce  peaks  in  Figures  15b  and  15c.  The  spectrum  with  a  distorted  interference 
peak  was  generated  by  a  z^,  in  which  the  J  interference  signals  on  the  right  hand  side  of  Eq. 
(30)  summed  destructively.  This  destructive  summation  was  also  the  reason  for  a  significantly 
stronger  signal  peak  in  the  spectrum,  because  the  denominator  on  the  right  hand  side  of  Eq. 
(27)  was  smaller. 

Figures  16a  to  16c  contain  the  LSAR  output  candidate  direction  estimates  calculated 
with  the  same  snapshots  used  in  Figures  14a  to  14c,  respectively.  To  avoid  overcrowding  in 
the  figures,  the  weakest  of  the  five  candidate  signal  direction  estimates  was  deleted  from 
every  snapshot.  As  SINR  increases  from  -50  dB,  Figure  16a,  calculated  with  C,  has  signs  of  a 
target  track  at  SINR=-35  dB  and  a  clearly  visible  track  for  SINR>-33  dB.  Figure  16b, 
calculated  with  CP,  has  signs  of  a  target  track  at  SINR=-18  dB  and  a  clearly  visible  track  for 
SINR>-15  dB.  It  has  false  tracks  at  0=-55°,  -7°  and  0°.  Figure  16c,  calculated  with  CGp,  has 
signs  of  a  target  track  at  SINR=-30  dB  and  a  clearly  visible  track  for  SINR>-27  dB.  It  has  a 
false  track  at  0=0°. 

The  above  results  show  that  the  use  of  C  instead  of  the  phase-only  Cp  could  reduce 
the  SINR  target  detection  threshold  by  approximately  17  dB,  when  an  extremely  strong 
interference  is  present.  Based  on  this  reduction,  the  use  of  C  could  increase  the  maximum 
nighttime  target  detection  range  of  the  HFSWR  by  over  50%.  The  results  also  show  that  the 
use  of  C  instead  of  the  gain-and-phase  Cgp  could  reduce  the  SINR  target  detection  threshold 
by  5  to  8  dB.  This  reduction  agrees  with  the  results  in  [15-18],  where  the  reductions  are 
between  4  and  8  dB. 


DFT  results 

The  DFT  estimates  of  signal  directions  and  powers  in  a  calibrated  array  snapshot  ycm 
were  calculated  with  the  following  procedure. 

Step  1. 

Calculate  a  new  snapshot  ycm  as 
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(32) 


yCm  = 


ymVN 

lyj 


Step  2. 

Specify  spectral  azimuth  direction  0. 

Calculate  the  array  steering  vector  a  for  this  direction  with  Eq.  (16),  using  0n,=0. 


Step  3. 

Calculate  the  complex  t)FT  amplitude  as 

z  =  iSk.  (33) 

N 


Step  4. 

Calculate  the  spectral  power  as 

P(0)  =  |zf  •  (34> 


Step  5. 

Repeat  Steps  2  to  4  with  many  azimuth  directions  equally  spaced  in  the  range  -90°  to 

90°  to  produce  a  DFT  spectrum. 

Identify  the  peak  positions  in  the  spectrum  as  candidate  signal  direction  estimates. 

Identify  the  peak  powers  as  signal  power  estimates. 

In  Step  1  the  average  power  at  the  array  elements  in  y  m  is  unity,  because  |y^jf=N  and  N  is  the 
number  of  array  elements.  In  Step  2,  the  division  by  N  ensures  that  z  is  equal  to  the  complex 
signal  amplitude  in  yCm,  if  only  one  signal  source  is  present  and  0  is  also  the  source  direction. 

Figures  17a,  17b  and  17c  show  the  DFT  spectra  calculated  with  the  synthetic 
snapshots  and  calibration  matrices  C,  CP  and  Cgp,  respectively.  The  target  and  interference 
peaks  are  significantly  broader  than  those  in  Figures  14a  to  14c.  There  are  many  false  peaks 
and  the  target  peak,  located  at  0=&r=-3O°,  is  weaker  than  the  nearest  false  peak  unless  SINR>- 
25  dB.  The  major  differences  among  the  figures  are  the  locations  and  the  power  levels  of  the 
false  peaks. 

Figure  18  shows  the  DFT  spectra  calculated  with  the  SINR=-50  dB  synthetic 
snapshot.  The  spectrum  calculated  with  C  (black)  has  eight  sidelobes  in  the  range  -75°  to  55°. 
It  looks  like  the  ideal  DFT  spectrum  of  a  pure  signal  with  source  direction  close  to  40°.  The 
spectrum  calculated  with  Cgp  (blue)  has  sidelobes  in  approximately  the  same  positions. 
However,  the  peak  powers  differ  noticeably  from  those  of  the  ideal  spectrum.  The  spectrum 
calculated  with  CP  (red)  only  has  six  sidelobes  in  the  range  -75°  to  55°.  The  sidelobe  positions 
and  powers  differ  noticeably  from  those  of  an  ideal  DFT  spectrum. 

Figures  19a  to  19c  contain  the  four  strongest  signal  direction  estimates  in  each  DFT 
spectrum  in  Figures  17a  to  17c,  respectively.  These  figures  do  not  have  a  firmly  established 
target  at  0=-3O°  unless  SINR>-25  dB.  There  are  at  least  two  false  targets  present.  The  major 


DREOTR  2001-103 


27 


differences  among  the  figures  are  the  scattering  of  the  target  direction  estimates  and  the 
locations  of  the  false  targets. 

The  above  results  show  that  it  is  difficult  use  the  DFT  to  detect  a  target  more  than  25 
dB  below  a  strong  and  highly  direction  interference  and  that  this  difficulty  cannot  be 
mitigated  with  an  accurate  receive  antenna  calibration.  The  results  are  not  unexpected, 
because  the  aperture  of  the  HFSWR  is  less  than  6  wavelengths  at  3.2  MHz,  the  DFT  has  eight 
sidelobes  in  the  range  -75°  to  55°,  and  the  peak  value  of  the  weakest  sidelobe  is  -24.0  dB 
below  the  main  peak. 
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Figure  6:  Leakage  power  calculated  with  C  (black),  phase-only  matrix  Cp  (red),  and  gain-and-phase 

matrix  Cgp  (blue). 
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Source  Azimuth  0m  (deg) 


Spectral  Azimuth  0  (deg) 


Figure  8a:  Directional  spectra  of  array  snapshots  in  preprocessed  calibration  data  file.  Matrix  C  was 

used  to  calibrate  the  snapshots. 
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Spectral  Azimuth  0  (deg) 


Figure  8c:  Directional  spectra  of  array  snapshots  in  preprocessed  calibration  data  file.  Gain-and- 
phase  matrix  CGp  was  used  to  calibrate  the  snapshots. 
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Figure  9a:  Directional  spectra  of  six  array  snapshots  with  source  azimuth  0m=4O°.  Matrix  C  was  used 

to  calibrate  the  snapshots. 
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Spectral  Azimuth  0  (deg) 


Figure  9c:  Directional  spectra  of  six  array  snapshots  with  source  azimuth  6m=40°.  Gain-and-phase 
matrix  CGp  was  used  to  calibrate  the  snapshots. 
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Figure  11a:  Directional  spectra  of  two  equal-strength  intersecting  targets.  Matrix  C  was  used  to 

calibrate  the  snapshots. 
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Spectral  Azimuth  0  (deg) 

Figure  11b:  Directional  spectra  of  two  equal-strength  intersecting  targets.  Phase-only  matrix  Cp 

used  to  calibrate  the  snapshots. 
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Figure  12b:  Some  directional  spectra  of  synthetic  snapshots  with  azimuth  0m=-3O°.  Phase-only  matrix 

Cp  was  used  to  calibrate  the  snapshots. 
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Spectral  Azimuth  0  (deg) 


Figure  12c:  Some  directional  spectra  of  synthetic  snapshots  with  azimuth  em=-30°.  Gain-and-phase 
matrix  CGP  was  used  to  calibrate  the  snapshots. 
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Figure  13a:  Candidate  direction  estimates  calculated  with  the  synthetic  snapshots  used  to  produce 
Figure  11a.  Matrix  C  was  used  to  calibrate  the  snapshots. 
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Source  Azimuth  0m  ( deg  ) 


Figure  13c:  Candidate  direction  estimates  calculated  with  the  synthetic  snapshots  used  to  produce 
Figure  11c.  Gain-and-phase  matrix  Cgp  was  used  to  calibrate  the  snapshots. 
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SINR  (dB) 


Spectral  Azimuth  0  (deg) 


Figure  14a:  Directional  spectra  of  a  target  in  strong  interference.  Matrix  C  was  used  to  calibrate  the 

snapshots. 
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Figure  14b:  Directional  spectra  of  a  target  in  strong  interference.  Phase-only  matrix  Cp  was  used  to 

calibrate  the  snapshots. 
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Figure  14c:  Directional  spectra  of  a  target  in  strong  interference.  Gain-and-phase  matrix  Cgp  was 

used  to  calibrate  the  snapshots. 
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Figure  15a:  Some  directional  spectra  of  synthetic  snapshots  with  SINR=-40  dB.  Matrix  < 

calibrate  the  snapshots. 
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Spectral  Azimuth  0  (deg) 

Figure  15b:  Some  directional  spectra  of  synthetic  snapshots  with  SINR=-23  dB.  Phase-only  matrix 

Cp  was  used  to  calibrate  the  snapshots. 
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Figure  15c:  Some  directional  spectra  of  synthetic  snapshots  with  SINR=-30  dB.  Gain-and-phase 
matrix  CGp  was  used  to  calibrate  the  snapshots. 
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Spectral  Azimuth  0  (deg) 


Figure  17a:  DFT  directional  spectra  of  array  snapshots  used  to  calculate  the  Capon  spectra  in  Figure 
14a.  The  snapshots  have  a  target  in  strong  interference.  Matrix  C  was  used  to  calibrate  the 

snapshots. 
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Spectral  Azimuth  0  (deg) 


Figure  17b:  DFT  directional  spectra  of  array  snapshots  used  to  calculate  the  Capon  spectra  in  Figure 
14b.  The  snapshots  have  a  target  in  strong  interference.  Phase-only  matrix  Cp  was  used  to  calibrate 

the  snapshots. 
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Figure  19a:  DFT  candidate  direction  estimates  calculated  with  matrix  C  and  the  synthetic  snapshots. 
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Azimuth  Estimate  0  (deg) 


Figure  19b:  DFT  candidate  direction  estimates  calculated  with  phase-only  matrix  CP  and  the 

synthetic  snapshots. 
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Figure  19c:  DFT  candidate  direction  estimates  calculated  with  gain-and-phase  matrix  Cgp  and  the 

synthetic  snapshots. 
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5.  Conclusion 


A  description  and  a  performance  evaluation  of  the  matrix  calibration  method  developed 
for  the  receive  antenna  array  of  the  HFSWR  at  Cape  Race,  Newfoundland,  have  been  presented. 
The  performance  evaluation  compared  the  matrix  method  with  two  other  methods:  (a)  the 
single-set  phase-only  calibration  method  being  used  and  (b)  an  assumed  single-set  gain-and- 
phase  method  that  equalizes  the  complex  element  output  amplitudes  of  signals  in  the  boresight 
direction.  Some  notable  observations  on  targets  within  55°  of  the  boresight  direction  are  listed 
below. 

(a)  The  leakage  of  signal  power  outside  the  source  direction  were  0.04%  and  1.29%  in  data 
sections  with  high  and  low  signal-to-interference-plus-noise  ratios  (SINRs),  respectively,  in 
calculations  with  the  matrix  calibration  method.  They  were  between  5.20  and  1 1.7%  in  the 
single-set  phase-only  method,  and  between  0.21  and  5.73%  in  the  single-set  gain-and- 
phase  method. 

(b)  The  errors  in  source  direction  estimates  were  usually  the  smallest  in  the  matrix  method  and 
the  largest  in  the  phase-only  method. 

(c)  The  errors  in  source  direction  estimates  had  a  component  linearly  dependent  on  the  source 
azimuth  in  the  phase-only  and  gain-and-phase  methods. 

(d)  The  signal  peaks  in  the  directional  spectra  of  the  array  snapshots  were  the  narrowest  in 
calculations  with  the  matrix  method.  They  were  significantly  broader  in  the  phase-only 
method.  Those  in  the  gain-and-phase  method  were  just  as  narrow  as  the  peaks  in  the  matrix 
method  for  signals  in  the  boresight  direction.  They  increased  to  about  half  the  widths  of  the 
phase-only  method  as  the  source  moved  towards  the  ±50°  directions. 

(e)  The  power  levels  of  false  peaks  were  less  than  -40  dB  of  the  true  signal  in  the  matrix 
method.  They  were  above  -30  dB  and  frequently  over  -20  dB  in  the  phase-only  method. 
Those  in  the  gain-and-phase  method  are  usually  between  -30  and  -20  dB. 

(f)  The  false  peaks  did  not  form  false  target  tracks  in  the  matrix  method.  They  frequently 
formed  false  tracks  in  the  phase-only  and  gain-and  phase  methods. 

(h)  Line  splitting  (two  direction  estimates  separated  by  less  than  one  beamwidth  in  high- 
resolution  direction-finding,  one  signal  source  actually  present)  was  very  common  in 
calculations  with  the  phase-only  method. 

(i)  The  SINR  threshold  in  detecting  targets  outside  the  region  of  extremely  strong  interference 
was  the  lowest  in  the  matrix  method.  It  was  approximately  17  dB  higher  in  the  phase-only 
method  and  5  to  8  dB  higher  in  the  gain-and-phase  method. 

Based  on  the  approximate  17  dB  difference  in  SINR  thresholds  for  target  detection,  the 
matrix  calibration  method  could  increase  the  present  maximum  nighttime  target  detection 
range  of  the  HFSWR  by  more  than  50%. 
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Appendix  A.  Calculation  of  signal  powers  and 
direction  estimates 


The  least  squares  autoregressive  method  [4,5]  was  combined  with  a  root-finding 
method  to  calculate  signal  direction  estimates  with  an  array  snapshot.  The  signal  power 
associated  with  each  direction  estimate  was  then  obtained  with  a  least  squares  fit  to  the 
snapshot.  In  the  procedure  below,  y  is  the  array  snapshot,  {yi,y2,...,yN}  are  the  components  of 
y,  N  is  the  number  of  array  elements,  K  is  the  order  of  linear  prediction  filter  in  the 
autoregressive  method  and  is  the  largest  integer  not  exceeding  N/3,  {gi,g2,—,gic}  are  the  filter 
coefficients,  {zi,z2,...,Zk}  are  the  poles  of  the  autoregressive  filter,  {0i,02,...,0k}  are  signal 
direction  estimates  calculated  with  the  poles,  {ai,a2,...,aic}  are  complex  signal  amplitude 
estimates,  and  {pi,p2,...,p]<}  are  signal  power  estimates.  There  are  six  steps  in  the  procedure. 
In  Step  1,  the  parameter  values  are  N=16  and  K=5.  In  Step  4,  a  is  a  vector  with  K  components 
and  is  the  least  square  solution  of  the  equation 

y  =  Aa  +  e,  (A.l) 

where  e  is  an  error  vector.  In  Step  5,  directions  with  |sin0iJ>l.O  are  non-physical  and  are 
ignored. 

Step  1. 

Solve  the  following  equation  for  the  coefficients  {ghg^—jgic}, 


u 

yK+i 

yK+i  yK  -  y2 

g. 

y*c+2 

... 

g2 

• 

yw-i  yN-2  —  yN-K 

* 

= 

yN 

y2  ys  •••  yic+i 

. 

y* 

... 

L  gK  J 

• 

L 


♦  * 

y  N-K+I  y  N-K+2 


Yn 


L  y*N-K  J 


Step  2. 


Find  the  roots  of  the  polynomial 
F(z)  =  zK-g,zK-,-g2zK-2-...-gK. 


Denote  them  by  {zi,z2,...,zK}. 


Step  3. 

Construct  a  matrix  of  array  steering  vectors  as 


(A.2) 


(A.3) 
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A  -  (ai,a2,...,aK) , 


(A.4) 


where  ak  is  an  array  steering  vector  given  by 
ak  =  (l,zkI,Zk,...,zkK*1)T, 

and  zk=zk/Jzk|- 


(A.5) 


Step  4. 

Calculate  a,  where  a=(ai,a2,...,aK)T,  as 

a  =  A#y .  (A.6) 

Here,  A#  is  the  Moore-Penrose  pseudo-inverse  of  A.  Use  the  singular  value 
decomposition  method  to  calculate  A#.  Equate  a  singular  value  to  zero  if  it  is  less  then 
8  times  the  largest  singular  value. 


Step  5. 

For  each  k,  calculate  sin0k  with  zk  and  the  relation 

2nd  sin0k  =  ang(zk) ,  (A.7) 

where  d  is  the  element  spacing  in  wavelengths  and  ang(Zfc)  is  the  phase  of  zk. 

If  Isine^l.O, 

Calculate  signal  direction  estimate  0k  with  Eq.  (A.7),  and 
Calculate  the  corresponding  signal  power  estimate  with  the  relation 


Pk  =  |akf . 


(A.8) 


Step  6. 

Arrange  the  direction  estimates  in  the  order  of  decreasing  signal  power  estimates. 


The  value  5=0.0001  was  used  in  this  report.  This  value  appeared  to  provide  a  good 
discrimination  between  signal  and  noise  subspaces. 
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Appendix  B.  Calculation  of  Capon  directional  spectra 


The  Capon  minimum  variance  method  [1-3]  is  combined  with  a  subarray  method 
for  uniform  linear  arrays  to  calculate  the  directional  spectrum  of  an  array  snapshot.  Two 
examples  are  then  given  to  clarify  the  procedure  and  the  definition 


_  Power  arriving  at  array  elements  from  spectral  direction  6  (B.l) 

*  Total  power  at  array  elements 

In  the  procedure  below,  y  is  the  array  snapshot,  {yi,y2,— O^n}  are  the  components  of  y,  N  is  the 
number  of  elements  in  the  full  array,  J  is  the  number  of  elements  in  the  subarrays,  K  is  the 

number  of  subarrays,  {zi,z2,...,zK}  are  the  subarray  snapshots,  {zi,z2,...,zK}  are  reversed- 

conjugated  subarray  snapshots,  e  is  a  noise  parameter,  a2  is  injected  isotropic  noise  power  at 

the  array  elements,  Ij  is  a  JxJ  identity  matrix,  a  is  the  full  array  steering  vector,  {ai,a2,...,aN} 
are  the  components  of  a,  and  b  is  the  subarray  steering  vector  constructed  with  the  first  J 

components  of  a,  i.e.. 


b  =  (ai,a2,...,aj)T  . 

The  parameter  values  are  N=16,  J=8,  K=N+1-J=9,  e=0.00001. 


(B.2) 


Step  1 . 

Construct  K  subarray  snapshots  as 


Zk  =  (yk,yk+i,-,yk-i+j)T , 


k=l,2,...,K .  (B.3) 


Step  2. 

Construct  the  reversed-conjugated  snapshots  as 

Zk  =  (y  k-i+j,.»»y  ku,y  k)T , 


k=l,2, 


„K. 


(B.4) 


Step  3. 

Construct  a  covariance  matrix  R  as 


R  =  (2K)-,|i(zkZkH  +  zkZkH}  . 


(B.5) 


Step  4. 

Calculate  o2  as 

o2  =  eTr(R) ,  C6-6) 
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where  Tr(R)  denotes  the  trace  of  R. 


Step  5. 

Construct  an  augmented  covariance  matrix  Ra  as 
Ra  =  R  +  c^Ij , 


(B.7) 


Step  6. 

Specify  azimuth  direction  0. 

Construct  the  array  steering  vector  a  for  the  full  array  with  Eq.  (16). 
Construct  the  subarray  steering  vector  b  with  (B.2). 

Step  7. 

Calculate  the  signal  power  from  direction  0  as 


P(0)  = 


1 

bHRAb 


1 

X  J-'Tr(RA)  ’ 


(B.8) 


The  product  J''Tr(RA)  is  the  total  power  at  the  array  elements  after  the  addition  of  the 
isotropic  noise  term  cPlj  to  R. 


Example  1.  One  signal,  no  noise,  no  interference 

The  signal  has  direction  0j,  amplitude  ai,  array  steering  vector  a,,  and  subarray 
steering  vector  bi.  The  snapshot  of  the  full  array  is  given  by 


y  =  a,a, , 

and  subarray  snapshot  zu  for  example,  is  given  by 


(B.9) 


Zi  =  aibi . 

One  can  now  verify  that 


(B.10) 


R  =  (otifb,bV , 


o2  =  eJ|aif , 


Ra  =  |aif(  bib?  +  eJIj) , 


Tr(RA)  =  W(l+eJ). 


(B.ll) 
(B.12) 
(B.  13) 
(B.14) 


DREOTR  2001-103 


71 


Using  the  matrix  inversion  lemma  and  noting  that  bVbi=J, 


R-i=  (l+e)J~'lj  -  J'2b|b? 

eO+eW 


(B.15) 


and 


prfn_  e(l+e)|aif 

(1+eH'W  ^  (l+eJJla.r 


(B.16) 


In  the  signal  direction,  where  0=8i  and  b=bi,  this  equation  gives 


P(6,)  = 


fop+elctil2 
laip+eJIaif  ' 


(B.17) 


The  numerator  on  the  right  hand  side  is  the  sum  of  the  signal  power  and  the  isotropic  noise 
power  from  the  signal  direction.  The  denominator  is  the  sum  of  signal  power  and  the  sum  of 
isotropic  noise  power  over  all  directions. 

In  directions  such  that  b^bi=0,  the  signal  is  absent  and  Eq.  (B.16)  gives 


P(0)  = 


efrif 

|a,f+EJ|a,p  ‘ 


(B.18) 


Here,  the  numerator  is  the  noise  power  from  direction  0.  In  the  limit  e=0,  the  right  hand  sides 
of  Eqs.  (B.17)  and  (B.18)  become  1.0  and  0.0,  respectively. 


Example  2.  Multiple  signal  sources 

There  are  J  signal  sources  with  directions  {0i,02,...,0j},  amplitudes  {ai,a2,...,aj}, 
array  steering  vectors  {ai,a2,...,aj},  and  subarray  steering  vectors  {bi,l>2,...,bj}.  To  simplify  the 
discussions,  two  assumptions  are  made.  First,  the  source  directions  are  mutually  orthogonal, 
so  that 


b^bq  =  J6pq  ,  (B.19) 

where  8pq  is  the  Kronecker  delta  function,  i.e.,  6pq=l  if  p=q  and  6pq=0  if  p*q.  Second,  K  is  so 
large  that  R  can  replaced  by  the  expected  value. 

In  this  example, 

y  =  ccia,  +  a2a2  + ...  +  ajaj ,  (B.20) 
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and  %\,  for  example,  is  given  by 


z,  =  aibi  +  a2b2  + ...  +  ajbj .  (B.21) 

Equating  R  to  the  expected  value  gives 

R  =  |otifb,b?  +  |a2pb2by  + ...  +  (ajfbjt^ ,  (B.22) 

0%  =  eJPsIj 

=  ePs(b,bV  +  b2b5+...  +  bJ^).  (B.23) 

where 

Ps  =  |aif  +  |a2f  + ...  +  lotjp  (B.24) 

is  the  sum  of  the  signal  powers  at  individual  array  elements.  Eq.  (24)  uses  the  property  that 

(b],b2,..,bj)**(bi,b2,..,bj)  =  JIj .  (B.25) 

This  leads  to 

JIj  =  (bi,b2,..,bj)(bi,b2,..,bj)H 

=  b,bV  +  b2l£  +  ...+bJltf .  (B.26) 

Substitutions  of  Eq.  (B.22)  and  (B.23)  into  (B.7)  and  using  Eq.  (B.19), 

Ra  =  flaif+ePs)b,bV  +  da2f+ePs)b2b?  + ...  +  fla,f+ePs)bj$ ,  (B.27) 

Ra1  =  flocf+ePsy^bibV  +  (Kxzf+ePsr'J’^by  + ...  +  flajf+ePsX^bjtf  ,  (B.28) 

Tr(RA)=  (l+eJ)JPs ,  (B.29) 

and 

P(ft)=  |(xji2+ePs  .  (B.30) 

1  S)  Ps+£JPs 

The  numerator  on  the  right  hand  side  is  the  power  at  an  array  element  due  to  the  signal  source 

at  direction  0j  plus  the  injected  isotropic  noise  power  from  this  direction.  The  denominator  is 
the  sum  of  signal  powers  from  all  sources  plus  the  sum  of  injected  noise  powers  over  all 

directions.  In  the  limit  e=0,  Eq.  (B.29)  gives  P(0j)=|(Xj p/Ps* 
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